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Abstract 

Finite time blow-up is shown to occur for solutions to a one-dimensional quasilinear 
parabolic-parabolic chemotaxis system as soon as the mean value of the initial condition 
exceeds some threshold value. The proof combines a novel identity of virial type with 
the boundedness from below of the Liapunov functional associated to the system, the 
latter being peculiar to the one-dimensional setting. 



1 Introduction 

We study the possible occurrence of blow-up in finite time for solutions to a one-dimensional 
parabolic system modeling chemotaxis [15]. More precisely, we consider the Keller-Segel 
chemotaxis model with nonlinear diffusion which describes the space and time evolution of 
a population of cells moving under the combined effects of diffusion (random motion) and a 
directed motion in the direction of high gradients of a chemical substance (chemoattractant) 
secreted by themselves. If -u > and v denote the density of cells and the (rescaled) con- 
centration of chemoattractant, respectively, the Keller-Segel model with nonlinear diffusion 
reads 

(1) dtu = div {a{u) Vu — uVv) in (0, oo) x f2 , 

(2) edtV = DAv--fv + u-M in (0, oo) x , 

(3) a{u) dyU = dyV = on (0, oo) x dVl , 

(4) (m,w)(0) = (uo,t;o) in fi. 
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In general, Q is an open bounded subset of M^, > 1, with smooth boundary dQ, a is 
a smooth non-negative function, and the parameters e, D, 7, and M are non-negative real 
numbers with D > and M > 0. In addition, the initial data Uq and vq satisfy 



(5) 



Mo > 0, 



/ 

Jn 



Uq{x) dx = \VL\ M , and 



/ 

Jn 



Vo{x) dx = . 



The constraints ([5]) ensure in particular that a solution {u, v) to satisfies (at least 

formally) the same properties for positive times, that is. 



The main feature of ([T]) is that it involves a competition between the diffusive term 
div {a{u) Vu) (spreading the population of cells) and the chemotactic drift term — div [u Vf ) 
(concentrating the population of cells) that may lead to the blow-up in finite time of the 
solution to The possible occurrence of such a singular phenomenon is actually an 

important mathematical issue in the study of which is also relevant from a biological 

point of view: indeed, it corresponds to the experimentally observed concentration of cells 
in a narrow region of the space which is a preamble to a change of state of the cells. From a 
mathematical point of view, the blow-up issue has been the subject of several studies in the 
last twenty years, see the survey [12] and the references therein. 

Still, it is far from being fully understood, in particular when e > (the so-called parabolic- 
parabolic Keller-Segel model). In that case, the only finite time blow-up result available 
seems to be that of Herrero & Velazquez who showed in [9l [10] that, when 1] is a ball in 
M^, D = 1, and a = 1, there are M > Sir and radially symmetric solutions {u,v) to 
which blow up in finite time. These solutions are constructed as small perturbations of time 
rescaled stationary solutions to and a similar result is also true when £ = [S]. The 

result in [10] actually goes far beyond the mere occurrence of blow-up in finite time as the 
shape of the blow-up profile is also identified. Recall that the condition M > Svr is necessary 
for the finite blow-up to take place: indeed, it is shown in [19] that, if f2 is a ball in M^, 
D = 1, and a = 1, radially symmetric solutions to ([I])-(11D are global as soon as M < Svr. 
We refer to [3, [IS] for additional global existence results when f2 is a bounded domain in 
M^, e > 0, and a = 1. In [TT] the existence of unbounded solutions is shown for e > and 
a = 1, but it is not known whether the blow-up takes place in finite or infinite time. The 
same approach is employed in [13] to obtain unbounded solutions to quasilinear Keller-Segel 
systems, still without knowing whether the blow-up time is finite or infinite. The finite 
time blow-up result proved in this paper (Theorem [T]) is thus the first one of this kind for 
quasilinear parabolic-parabolic Keller-Segel systems. 

In contrast, for the parabolic-elliptic Keller-Segel system corresponding to e = 0, several 
finite time blow-up results are available. There is thus a discrepancy between the two cases 
e > and e = which may be explained as follows. On the one hand, as observed in [13] 



(6) 




and 
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when e = 0, f2 is a ball of M^, a = 1, and uq is radially symmetric, it is possible to reduce 
([T])-(jl]) to a single parabolic equation for the cumulative distribution function 



Finite time blow-up is then shown with the comparison principle by constructing appropriate 
subsolutions. This approach was extended to nonlinear diffusions (non-constant a) and 
arbitrary space dimension > 1 in [6] . On the other hand, it has been noticed in [21 [16] 
that, still for a = 1, the moment of u defined by 



satisfies a differential inequality which cannot hold true for all times for a suitably chosen 
value of > 0, for it would imply that u reaches negative values in finite time in contradiction 
with In contrast to the previous approach, this is an obstructive method which provides 
no information on the blow-up profile and is somehow reminiscent of the celebrated virial 
identity available for the nonlinear Schrodinger equation (see, e.g., [H Section 6.5] and the 
references therein). Nevertheless, it applies to more general sets Q [TTl [181 [20] . We recently 
develop further this technique to establish finite time blow-up of radially symmetric solutions 
to ([ID-(|1D with e = in a ball of M^, N >2, when the diffusion is nonlinear [S], the main idea 
being to replace the moments by nonlinear functions of the cumulative distribution function 
U. For a related model in M.^ with nonlinear diffusion a{u) = m m > 1, finite time 

blow-up results were recently established in [3l [21] by looking at the evolution of the second 
moment M2. 

Coming back to the parabolic-parabolic Keller-Segel system (e > 0), it seems un- 

likely that the first approach described above (reduction to a single equation) could work 
and the purpose of this paper is to show that finite time blow-up results can be established 
by the second approach in the one- dimensional case (A^ = !)• More precisely, we consider 
the initial-boundary value problem 

(7) dtu = dx {a{u) d^u — u d^v) in (0, 00) x (0, 1) , 





k e (0, 00) , 



(9) 
(10) 
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e dtv 



Ddlv -'J v + u- M 
a{u) dxU = dxV = 
iu,v){0) = {Uo,Vo) 



on 



m 



m 



(0,cx)) X (0,1), 
(0,cx)) X {0,1}, 
(0,1), 



and assume that 



(11) 



£>0, D>0, 7>0 



M > 0, 



and the initial data (uq^Vq) G W^''^{0, 1;M' 



.^) satisfy 



(12) 




and 
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We further assume that a G C^(M) and that there are p E (1,2], and ci > such that 

(13) < a(r) < ci (1 + r)-P for r>0. 
Our main result then reads as follows. 

Theorem 1 Assume that the parameters e, D, 7, M, the initial data {uq,vo), and the 
function a fulfil the conditions f|TT]) . f|T2l) . and f|T3|l . respectively. Then there is a unique 
classical maximal solution 

{u,v)eCi[0,T^) X [0,l];M2)nCi'2((0,Tj x [0,1];^^) 
to (l71)- (fT0l) with maximal existence time G (0, 00]. It also satisfies 

(14) M(t,a;)>0, I u{t,x)dx = M, and f v{t,x)dx = 

Jo Jo 

for {t,x) G [0,Tm) X [0,1]. Introducing 

F{zi, Z2) := ci {l + M) + ^ + zi + M Z2 + ^^-^ zl , 

(15) P,(zi, ^2, ^3) := y + ^ + + 1^ ^3 j F{z,, Z2) 



+ TTTTrn F(zi,Z2) 



(p- i)Mp-i ' ' ^ g(g + 1) 

and 

1 / /"x \ g 

^^0(2/) dy) dx, 



/or {zi,Z2,Z3) G [0,00)^ anc? g > 2, ti;e /iawe < 00 as soon as Vq{mg{0), \\vQ\\Hi,eM) < 
for some finite g G (2, 2/(2 — p)]. In particular, if Uq is such that 

(16) (mg(0),0,0) < for some finite g G (2, 2/(2 - p)] , 

i/iere zs ^9 > such that eM G (O,-;?) and ||fo||//i < im;)/^/ that Vq (mg(0), ||i'o||_f/i, £^) < 
and thus Tm. < 00. 

There are functions uq satisfying ( fT2|) and ( fT6l) if M is sufficiently large. Indeed, observe 
that 



P,(0,0,0) = (1 + ^) (^ci(l + M) + ^^ 



(p-l)MP-i V ' 2D J q{q + K 
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is negative for sufficiently large M as g > 2. Given such an M > and choosing the 
function uo(x) = 2M max{x + 5- l,0}/5^, x E (0, 1), we have mg(0) = (2M)«(5/(2g + 1) 
and (mg(0), 0, 0) < for 5 > small enough. In fact, if Uq fulffis f|T6|) . then the same 
computation as the one leading to Theorem [T] shows that the corresponding solution to the 
parabolic-elliptic Keller-Segel system {e = 0) blows up in a finite time and the last assertion 
of Theorem [1] states that this property remains true for the parabolic-parabolic Keller-Segel 
system (e > 0) provided e and vq are small, that is, in a kind of neighbourhood of the 
parabolic-elliptic case. 

Remark 2 The growth condition required on a in f|T3|) is seemingly optimal: indeed, it is 
proved in that Tm. = oo if a{r) > cq (1 + r)~^ for some p < 1 and e = 0, and the proof is 
likely to extend to the case e > 0. 

The proof of Theorem [1] relies on two properties of the Keller-Segel system (17j)- (fT0l) : ffist, 
there is a Liapunov functional [7] which is bounded from below in the one-dimensional case 
[6] and which provides information on the time derivative of v. This will be the content of 
Section [2] where we also sketch the proof of the local well-posedness of (l7l)- fllUI) . We next 
derive an identity of virial type for the L'^-norm of the indefinite integral of u in Section [3] 
which involves in particular the time derivative of v. The information obtained on this 
quantity in the previous section then allow us to derive a differential inequality for the L"^- 
norm of the indefinite integral of u for a suitable value of q which cannot be satisfied for all 
times if the parameters e, D, 7, M, and the initial data [uq, vq) are suitably chosen. 



2 Well-posedness and Liapunov functional 

In this section, we estabhsh the local well-posedness of fl7|)- (fT0l) in W^''^{0, 1;M?) and recall 
the availability of a Liapunov functional for this system [7]. To this end, we assume that 

(17) < a e C\R) 
and define b G C^((0, 00)) by 

(18) 6(1) = 6'(1) := and b"{r):=^ for r>0. 

r 

Proposition 3 Assume that the parameters e, D, 7, M, and the function a fulfil ( fTTl) and 
(ITTl) . respectively. Given initial data (uo,fo) G Vr^'^(0, 1; M^) satisfying ( |T2l) . there is a 
unique classical maximal solution 

iu,v)eCi[0,T^) X [0,l];M2)nCi'2((0,Tj x [0,1]; M^) 

to (171)- (fT0l) with maximal existence timcTm G (0, 00] and {u,v) satisfies f|T4|) fort G [0,Tm). 
In addition, ifTm < 00, we have 

(19) lim (||M(t)|U + ||t;(t)||oo) = 00. 

t — >Tm 
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Proof. We define a G C'^{W^]M2{ 

d{y) = (5'"'"(2/))i<^,,<2 := I I for y = (^1,^2) G 

and introduce tlie operators 

A{y)z := -d^ {d{y) d^z) , 
{B{y)z{f)),B{y)z{l)) := {-~a{y) d^z{f))ra{y) d^z{l)) 




-1 yi + y2- M 



and tlie function 

m ■■- 



witli z = {zi, Z2). Witli tliis notation, an abstract formulation of fl71)- f|T0l) reads 

dtz + A{z)z = f{z), 
B{z)z = 0, 

z{0) = {vo,Uo), 

with z = {v, u) and we aim at applying the theory developed in [1]. Owing to (fTT]) and ([1] 
d{y) is a positive definite matrix for all ?/ G and we infer from [1], Section 4] that the 
boundary-value operator {A., B) is normally elliptic. It then follows from [H Theorems 14.4 
& 14.6] that (171)-( |T0l) has a unique classical maximal solution 

{v,u) :=zgC([0,T„) X [0,l];M2)nCi'2((0,T„) x [0,1]; M^) 

for some G (0, 00]. In addition, a^'^(?/i, 0) = for i/i G M and we deduce from [H Theo- 
rem 15.1] that u(t,x) > for {t,x) G [0, T^) x [0, 1]. The property (fT4|) then readily follows 
from (17|)- (ITU]) and (IT^ by integration. As for the last statement (HM . it is a consequence of 
the lower triangular structure of the matrix a and [H Theorem 15.5]. □ 

We next proceed as in |^ to check the availability of a Liapunov functional for fl7|l- f|T0|) . 



Lemma 4 Assume that the parameters e, D, 7, M, and the function a fulfil f[TT|) and 
f|T7|) . respectively. Given initial data (uo,t'o) G W^''^{0,1]M.'^) satisfying f|T2|) and such that 
b{uQ) G -/^""^(0, 1), the corresponding classical solution {u,v) to fl7|)- ffTI]|) satisfies 

(20) / < L(mo,^^o) M tG[0,T^), 

Jo 

where 

(21) L{u,v):= [ ( b{u) - uv + ^ Id^vl"^ + \v\A dx . 
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Proof. It follows from ©-Q that 



—L{u,v) = / {b'(u)—v) dtU dx + / (D dxV dxdtv + {•y v ~ u) dtv) dx 
dt Jo Jo 

1 

{h"{u) dxU — dxv) [a{u) dxU — u dxv) dx 





/ dtv {—D d^v + 7 y — u) dx 
Jo 



-1 2 
u \dx (b' (u) — v)\ dx — (M + e dtv) dtV dx 
'0 Jo 

(22) < -s\\dM\l 

the last inequality being a consequence of (fT4l) . Integrating the previous inequality with 
respect to time gives (l20l) . □ 

We next take advantage of the one-dimensional setting to show that L is bounded from 
below without prescribing growth conditions on a. This fact has already been observed in 
[6] and is peculiar to the one-dimensional case. Indeed, as shown in [TJ [11] , the occurrence 
of blow-up is closely related to the unboundeness of the Liapunov functional. 

Lemma 5 Assume that the parameters e, D, 7, M, and the function a fulfil f|TT]) and 
f fT7|) . respectively. Given initial data {uq,vo) G 1V^'^(0, 1; M^) satisfying f|T2|) and such that 
b{uo) G L^(0, 1), the corresponding classical solution {u,v) to (l71)-(TT0l) satisfies 

(23) L{u{t),vit))>-^ for t G [0, TJ . 

Proof. Owing to f|T4l) . the Poincare inequality ensures that ||t'(t)||oo < llf^i^f (t)||2 for t G 
[0, Tm) so that 

1 

uit) v{t) dx< \\v{t)\\^ \Ht)\\i < \\dxv{t)\\2 \\u{t)\\i. 



We use again (fT4l) as well as the non-negativity of h to conclude that 

L{u{t)Mt)) > § WdMml - M \\dxvm2 = y {\\dxvm2 - f )' - ^ 
for t G [0,Tm), from which ( l23i) readily follows. □ 
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3 Finite time blow-up 



As already mentioned, the main novelty in this paper is a new identity of virial type which is 
the cornerstone of the proof that blow-up takes place in finite time under suitable assump- 
tions. Specifically, we assume that the parameters £, 7, M, and the function a fulfil the 
conditions ffTT]) and ffT5]l . respectively. Recalling the definition (ITSl) of 6, we deduce from 
(DS]) that 

(24) h{r) < Cl (r Inr - r + 1) l[o,i](r) + ~ l[i,oo)(r) < Ci (1 + r) , r > . 
We also define 

a{s)ds, r>0, 

and infer from (fT3l) that A is well-defined and satisfies 

(26) < -Air) r < r^-^ , r > . 

p — \ 

Consider next initial data {uq.Vq) G iy^'^(0, 1; M^) satisfying (fT2|) . If (M,f) denotes the 
corresponding classical solution to (171)- (fT0l) given by Proposition [3l we define the cumulative 
distribution functions U and V by 



(27) U{t,x):= u{t,y)dy and V{t,x) := v{t,y) dy 

Jo Jo 

for {t,x) e [0,T^) X [0, 1]. It readily follows from and ([HD that ([/, \^) solves 

(28) dtU = d,A{u) - u d,v in (0, T„) x (0, 1) , 

(29) e dtV = D d^v - V + U - Mx in (0, T^) x (0, 1) , 

the function A being defined in (125!) . and 

(30) U{t,0) = M -U{t,l) =0 and ^(t, 0) = ^(t, 1) = , t G [0, T^) . 
Lemma 6 Introducing mq(t) := ||f/(i)||q/q' /or q > 2, we have 



dm M Mi+^ 

(31) ^ = -p - _^ + M"-^ A(^(t, 1)) - (g - 1) U^'~' u A{u) dx 

+ — / t/« 9^1; rfx - — / f/^~^ uV dx 

qD Jo D Jo 

forte [0,T^). 
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Proof. We infer from ([28]), ([29]), and ([30]) that 

dmg _ r „_i A,^-\x=i I X TTq-2 



1 

- — / n U''^^ {edtV + -fV -U + Mx) dx 
" Jo 

= M''-' A{u{t, 1)) - (g - 1) / U'^^^ u A{u) dx-^ [U" dtV]lZl 

Jo 

h [ * - i [ '^^"^^ [^"']"; 



qD 

qD ^ J-=o ^ D ' 



M^~^ A(u(t, 1)) - (g - 1) [ V^-^ u Aiu) dx + ^ I W dtv dx 

Jo qD Jo 



^ I U"-^ uV dx- -^^^^-^vt: + ^ , 



D J, ^ ' q{q + l)D + 'D 

which is the expected identity. □ 



At this point, we notice that the solution to the ordinary differential equation D X = 
MX — [M'^^^ / [q[q + 1))) (obtained by neglecting several terms in fl3ip ) is given by 

and thus vanishes at a finite time if X(0) < M'^/{q{q + 1)). If a similar argument could 
be used for ( I3T]) . we would obtain a positive time to such that mq(to) = which clearly 
contradicts the properties of Uito): indeed, by (1271) and (130]) . x i— U(to,x) is continuous 
with U{tQ, 1) = M. Consequently, the solution {u,v) to (17))- (fT0]) no longer exists at this time 
to and blow-up shall have occurred at an earlier time, thus establishing Theorem [H For 
this approach to work, we shall of course control the other terms on the right-hand side of 
( 13T]) which will in turn give rise to the blow-up criterion stated in Theorem [H The latter is 
actually a simple consequence of the following result: 

Theorem 7 Assume that the parameters e, D, 7, M , and the initial data (uo, vq) are such 
that 

( M^\ 
(32) E m,(0) + L{uo, vo) + 777T < 



2D J 

for some finite q E {2,2/(2 — p)], where 



DM" AqD J (p-l)MP-i g(g + 1) 

for 2 > 0. Then < 00. 
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Proof. The starting point of the proof being the identity fl3T]) . we first derive upper bounds 
for the terms on the right-hand side of fl3Tl) involving A, e, and 7. Thanks to fl26|) and the 
non-negativity of U, it follows from the Holder inequality that 



•1 

M"-^ A{u{t, 1)) - (g - 1) / f/^"^ u A{u) dx 
P - 1 io 



< 



Since g G (2, 2/(2 — p)], we may use the Jensen inequality and (fT4|) to conclude that 

(33) ikf-i Aiuit, 1)) - (g - 1) u Aiu) dx < ^^(g-^)^^' ]\^2-p ^(9-2)/^ _ 

io - 1) 

Next, to estimate the term involving 7, we adapt an argument from [16] and first claim that 

(34) Vit,x)>Vm{t,x):=^ ix'-x) + hit,x), (t, x) G [0, TJ x [0, 1] , 
where h denotes the unique solution to 

(35) e dth- D dlh + -f h = 0, (t, x) G (0, 00) x (0, 1) , 

(36) h{t,0) = h{t,l) =0, tG(0,cx)), 

(37) h{0, x) = min |y(0,x) + ^ (x-x3),o| < 0, xG(0,1). 

Indeed, Kn < ^ on [0,Tm) x {0, 1} and {0} x [0, 1], and it follows from the non-negativity 
of U and the negativity of h that 

edtV^,-DdlV,n + lV,n = edth-Mx-D dlh + ^ {x^ -x)+-i h 

< -Mx <U - Mx = e dtV - D dlV + V . 

The comparison principle then implies flMl) . We next infer from 0341) and the non-negativity 
of u and U that 

_2 f f/9-1 uV dx < ! f/""^ u Vr, dx 



D Jo - D JO 



qD qD Jo 



7 / M 
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We next note that dxh also solves fl35l) with homogeneous Neumann boundary conditions, 
the latter property being a consequence of fl35l) and fl36l) . Since 



|a^/i(o,x)| < 



II II M 

< IK'o||oo + 



the comparison principle and the non-negativity of 7 warrant that ||52;/i(t)||oo < ll'^olloo + 
{M/3D) for t > 0. Consequently, recalling the Sobolev embedding ||fo||oo < IK'olli^i) we end 
up with 



(38) 



1_ 
D 







uV dx<^ 



D 

M 



\vo\\m "1, 



We finally infer from (fT4l) . (127|) . (l30l) . and the Holder inequality that 



(39) — / U" dtv dx < 

^ ^ qD Jo ' - qD 



Jo q 



It now follows from ([3l]), dM]), ([38]), and ([39]) that 



dm„ M 
* < — 
- D 



dt 



ci(g- l)g(''-2)/'?L) 
{p - l)MP"i 



{g-2)/g 



g(g + l) 



^^'^^ 1/2 n o „ 



< 



M 



Owing to (1121) and we have b{uo) G L^(0, 1) and it follows from 

above inequality that 



and the 



d 
dt 



rria + L(u, v) + 



M2\ 



M 



2D J 



I O II 1/2 



Using now the monotonicity of E and ([23]) . we end up with 

dt + + 2D j - D i"^^ + + 2D 



Assume now for contradiction that = 00. The previous inequality and fl32]) then warrant 
that there is a time to > such that mqito) + L(M(to), t'(to)) + (M^/2D) = and hence 
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TTT'qito) = by fl23|) . This in turn implies that U{to,x) = for all x G [0, 1] and contradicts 
fl30|) . Consequently, < oo. □ 

The remaining step towards Theorem [1] is to use the properties of a to simplify the condition 
fl32|) derived in Theorem [71 

Proof of Theorem [T]. It follows from ([T^ . ([21]), and the Sobolev embedding ||wo||oo < 
1 1 'Uo 1 1 HI that 

N /"V / X ^ ,o ,2 7 , ,2 n n \ , 

L(Mo,fo) + ^ < J I ci (1 + no) + — |d:r^^o| +^\vo\ +%|K'o||ool "^ + ^ 

< ci (l + M) + ^ + ^^^ ||t;o||^i + M llnoll^i 
= F(mg(0), llnollHi) - "^9(0) 5 
the function F being defined in Theorem [TJ Therefore, 

E f m,(0) + L(no, ^o) + ^ j < o F) K(0), \\vo\\m) = V, K(0), ||i;o||^,i, eM) , 

and the condition Vq (mq(0), ||fo||/i-i, eM) < clearly implies ([321) and hence Tm < 00. □ 
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Institut de Mathematiques de Toulouse, Universite Paul Sabatier. T. Cieslak would like to 
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